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Abstract. Consider the generalized Kuramoto-Sivashinsky (gKS) equation. It is a model pro- 
totype for a wide variety of physical systems, from flame-front propagation, and more general front 
propagation in reaction-diffusion systems, to interface motion of viscous film flows. Our aim is to 
develop a systematic and rigorous low-dimensional representation of the gKS equation. For this 
purpose, we approximate it by a renormalization group equation which is qualitatively characterized 
by rigorous error bounds. This formulation allows for a new stochastic mode reduction guaranteeing 
optimality in the sense of maximal information entropy. Herewith, noise is systematically added to 
the reduced gKS equation and gives a rigorous and analytical explanation for its origin. 

These new results would allow to reliably perform low-dimensional numerical computations by 
accounting for the neglected degrees of freedom in a systematic way. Moreover, the presented reduc- 
tion strategy might also be useful in other applications where classical mode reduction approaches 
fail. 

Key words. Generalized Kuramoto-Sivashinsky equation, renormalization group method, 
stochastic mode reduction 

Subject classifications. 35Q53, 37F25, 81T15, 81T17 

1. Introduction We consider abstract evolution equations of the form, 

u(x,0) = u o {x), 

where A denotes a general linear operator and B represents a nonlinear term of Burg- 
ers' type. Well known equations in this class include, e.g. the viscous Burgers equa- 
tion, the Korteweg-de Vries equation, and the Benney-Lin equation. We start by per- 
forming a formal renormalization group (RG) approach for the general form in 
We subsequently focus on a rigorous low-dimensional reduction of the generalized 
Kuramoto-Sivashinsky (gKS) equation: 

d t u + \uu x + Ku xx + Su xxx + isu xxxx = Q 1 m.V a x]0,T[, 

u(x,0)=g(x), in? Q , 

where V a :=] — air, air{ is a periodic domain with a := -k- for an arbitrary period L > 
while the solution u(x,t) : "P Q x]0,T[— of (|1.2|) describes for example the height of 
a one-dimensional surface above a substrate point x at time t. We also take g(x) € 
H g (V a ) for g>4, a periodic initial condition, i.e., 

g(x+L)=g(x). 
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Renormalization of the gKS equation 



The gKS equation is of the type (jl.ip with, 



B(m) := B(u,u) := Xuu 



and A:=(/sfi£+<yag + i/S£). 



(1.3) 



A rigorous dimension reduction of the gKS equation is of special interest because it 
is not a Hamiltonian system and does not have an intrinsic invariant measure. This 
makes direct application of stochastic mode reduction strategies difficult, see [33J for 
instance. 

It is noteworthy that the gKS equation retains the fundamental elements of any 
nonlinear process that involves wave evolution: the simplest possible nonlinearity 
uu x , instability and energy production u xx , stability and energy dissipation u xxxx and 
dispersion u xxx . We notice that the nonlinearity arises effectively from the nonlinear 
correction to the phase speed, a nonlinear kinematic effect that captures how larger 
waves move faster than smaller ones. In the context of thin-film flows [6, 361(37], the 
terms uu Xl u xx , u xxx and u xxxx are due to the interfacial kinematics associated with 
mean flow, inertia, viscosity and surface tension, respectively, with the corresponding 
parameters, A, k, 8 and v all positive and measuring the relative importance of these 
effects. The "strength of the nonlinearity", A, in particular, is associated with the 
scaling for the velocity (and hence time). In addition, j v udx — UQ, a measure of 
the volume of the liquid, a conservation property for systems whose spatial average 
does not drift. A simplified form of (|1.2[) is obtained by appropriately rescaling u, x 
and t which is equivalent to setting A = k = v = 1 and keeping the same notation for 
dimensionless quantities [5] [35] |37] . 

As with many nonlinear time-dependent problems in science and engineering, 
equations of the form (jl.lj) are too complex to be fully resolved and the influence of 
neglected degrees of freedom is not clear a priori. This problem exists independently 
of spatial dimensions for (|1.1[) and hence for the gKS equation also. The reliable 
resolution of high dimensional problems is a well-known issue in computational science 
where one can numerically only deal with a finite number of degrees of freedom. 

Hence, there is a strong need for dimensionally reduced formulations, which in 
turn would allow for studies of long time behavior of physical systems. Modeling of the 
ocean-atmosphere, which mainly generates our weather, is one important example: 
One has a characteristic timescale of several years for the ocean in contrast to a 
couple of days governing atmospheric structures such as cyclones. As a consequence, 
a characteristic feature of many physical systems is the presence of fast and slow 
degrees of freedom. The relevant information of a system's long time behavior is 
often primarily contained in the slow modes. 

This immediately leads to the general, often ad hoc, approximation of decompos- 
ing the problem of interest into fast w and slow v modes. For equation (|1.1|) such 
a purely formal splitting, i.e., u ~ u f ~ — v e + w e , reads in standard notation applied in 
the literature as, 



where the small parameter < e <C 1 mediates the timescale separation. Mode re- 
duction strategies, such as "adiabatic elimination" [35], invariant manifolds [7], and 
optimal prediction [3J , are tools to eliminate the fast modes and derive "appropriate" 




f(v e ,w e ), 



(1.4) 
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equations for the slow modes only. We remark that especially for systems with spatio- 
temporal chaos (like the gKS equation [HIM]) such a reduction needs to be carefully 
performed in order to not lose the relevant dynamical characteristics of the full system. 
The strategy of defining an invariant manifold is almost classical by now. For exam- 
ple, in [8] the existence of an inertial manifold for the KS equation (obtained from 
(|1.2p with 5 = 0) is shown. An inertial manifold is a finite-dimensional, exponentially 
attracting, positively invariant Lipschitz manifold. The principle idea is to determine 
a map $ : V — > W such that we can rewrite equation in the low-dimensional form, 

d t v + PB{v + $>{v),v + ${v)) + APv = 0, (1.5) 

where P:H^V and Q:=(I-P):i->W are projections onto the orthogonal sub- 
spaces V and W such that H = V®W. A strategy to determine $ in general Galerkin 
spaces is, for example, suggested in [7] for the KS equation. The RG approach per- 
formed here can also be understood as a formal and feasible procedure to derive an 
asymptotic invariant manifold, see (|1.8I) and (ll.9p . 

Open questions and answers to the classical separation (jl.4|) : (i) Is the 
splitting |J.^[ ) and the approximation of u by u e valid and in which sense? This ques- 
tion is often not answered in the literature where from the outset a separation (|1.4D 
is assumed, see [TJ 0] for instance. These studies heuristically motivate a timescale 
mediation (or separation into slow and fast scales) of the form (|1.4j) . The present work 
aims to provide a rigorous foundation in Theorem 13.61 by the following estimate, 

\\u-vX* {Va )( T )<Ce* + exv{CT)(e l l* + e) . (1.6) 

If we suppose that u and u e satify a Gevrey regularity characterized by a parameter 
cr > 0, then we can improve (|1.6p in the following way, 

ll«-^l!^(^)( T )<^ 2 + cx p(CT)( el/4ex p(-^r/4)+ e ))- (1-7) 

It should be pointed out that these estimates also account for the reduction to the 
slow degrees of freedom v e too and not only for the error between u and u £ . 

(ii) How can we account for the fast degrees of freedom w t in an equation for the 
slow modes v e only? For this purpose we apply an abstract RG approach extended 
to general multiscale problems, see O [24l [25] . The RG method was first introduced 
in quantum field theory as a tool to perform scale transformations. The method 
then became popular with Wilson's work on the Kondo problem |39j . It can formally 
provide the separation (|1.4p . This means, we first obtain an approximation for v e of 
the form 



d t v e + A v v* + P N B{v e ,v t ) = -eG e (U{t),v e ) , (1.8) 

where eG t {U,v e ) is a perturbation "force" originating from the renormalization 
method and U = V + W = v e + W is a solution of the RG equations 

d t V + A v V + P N B{V,V) = 0, 

d t W+Q N B x {V,W) = Q, ( ' ' 

where V — PnU and W = (I — Pn)U=: QnU are projections onto the slow and fast 
manifolds, respectively. Since we can analytically solve for W, we end up with an 
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Renormalization of the gKS equation 



equation for the slow variable v e only. The above estimates (|1.6p and (|1.7[) then make 
the reduction (jl.8p rigorous. Moreover, equation (|1.9p ^ can be interpreted as the map 
$ e (?/) onto the asymptotic invariant manifold. 

It should also be pointed out that at this stage the RG approximation (|1.8I) alone 
is not satisfactory since the fast variable W contained in U is of infinite dimension and 
hence can not entirely be resolved numerically. We give an answer to this problem 
after the last question (iii) by the principle of maximum entropy. Moreover, question 
(ii) is of particular relevance here, since the fast modes prevent the existence of a 
canonically invariant measure. Such a measure is the key information for classical 
reduction methods such as Mori-Zwanzig and optimal prediction, see question iii). 

(iii) What kind of information do we need to carry over from the (infinite dimen- 
sional) fast degrees of freedom to the (finite dimensional) slow ones and how? To 
this end, we derive a stochastic evolution equation for the resolved (slow) variable 
by properly including necessary information from the unresolved (fast) variable by 
a maximum information entropy principle introduced in |141 1151 128) . This principle 
does not require statistical data to define all Fourier modes. It turns out that the 
asymptotic behavior in time of a weighted variance of the fast modes is sufficient. 
The necessity of such a strong assumption relies on the fact the gKS equation does 
not have an infinite-dimensional invariant measure and that we only account for spa- 
tial randomness. Via this entropy principle (Theorem 14. ip we then conclude that 
the Fourier modes of the fast variable W in fOJh are Gaussian distributed with zero 
mean. Hence, we rigorously obtain a noisy gKS equation by applying the random vari- 
able U = v e + W in the deterministic equation for the slow variable 1(1.80 ■ Herewith, 
our analysis explains how to rigorously add a random force to the gKS equation. 
Furthermore, our derivation further shows that the induced noise accounts for the 
unresolved degrees of freedom and hence becomes less important for an increasing 
number of grid points in computations. 

The approach proposed here provides an alternative to the Mori-Zwanzig for- 
malism (26j HOI Hi] which depends on Hamiltonian 26, 40] or extended Hamiltonian 
structures (41) . Mori-Zwanzig techniques and related optimal prediction methods [3] 
generally rely on a canonical probability distribution (invariant measure) which exists 
naturally for Hamiltonian systems. The canonical probability density for a Hamilto- 
nian H(u) is p(u) :— Z _1 exp(— f3H(u)), where (3 is the inverse temperature and Z a 
normalization constant referred to as the partition function. The Mori-Zwanzig for- 
malism then is based on a projection operator P that projects functions in L 2 onto a 
subspace that only depends on the resolved degrees of freedom. With respect to the 
canonical density p such a projection operator P can be defined by the conditional 
expectation 

[P/](,):=E[/l,]= ^(;7>y^ , (1.10) 
1 J p[v,w) aw 

where /6L 2 and v is the resolved and w the unresolved variable. The projection P 
and Dyson's formula for evolution operators then provide an equation for the resolved 
modes v only. Moreover, (|1.10|) is the conditional expectation of / given v and hence is 
the best least square approximation of / by a function of v. Therefore, the projection 
P guarantees optimality which is the key idea in the optimal prediction method. 
However, neither a Hamiltonian structure nor an invariant measure exists for the gKS 
equation. Therefore, it is not obvious how to derive standard optimality statements 
relying on a conditional probability argument (33) . In contrast to such a conditional 
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probability approach, we achieve optimality in the sense of maximum information 
entropy. However, we remark that one can also define other projections than (|1.10[) . 

The purpose of the present article is threefold: 1. To reliably perform a (stochas- 
tic) mode reduction for the full gKS equation in contrast to [33 where a truncated 
problem is studied. The principal idea is based on an abstract RG approach, as 
emphasized earlier. We derive error estimates (Theorem I3.6j) for this reduction and 
hence provide rigorous support for the heuristic motivation of a noisy, low dimensional 
approximation deducted in [4] by the standard RG method in physics |39) : 

2. To rigorously support Stinis' assumption of Gaussian distributed Fourier modes 
[33] . To this end, we derive a probability distribution (Theorem I4.ip for the fast 
modes by the principle of maximum information entropy. 

3. The findings in 1 and 2 form the bases for a new stochastic mode reduction strat- 
egy. We arc able to reduce the fast variable by an equation for the slow variable only. 
The information of the fast modes enters as a random variable W via a force term 
into the slow mode equations. We are not aware of any previous work that utilizes 
the RG method in the context of stochastic mode reduction. 

We introduce basic notation and well-known results in Section 11.21 A formal 
derivation of an RG equation for the gKS equation follows in Section [2j In Section 
[3] we obtain error estimates to rigorously verify the approximation derived in Section 
[2J In Section [D we reduce the fast modes by a mode reduction strategy based on 
the maximum information entropy principle. Finally, in Section [6] we close with 
conclusions and perspectives. 

1.1. The gKS equation The KS equation is a paradigmatic model for the 
study of low-dimensional spatio-temporal chaos or weak/dissipative turbulence as de- 
fined by Manneville |23j . This type of turbulence is often characterized by formation 
of clearly identifiable localized coherent structures in what appears to be a randomly 
disturbed system, as is e.g. the case with Rayleigh-Benard convection [3Tj . The KS 
equation was first proposed as a model for pattern formation in reaction-diffusion 
systems by Kuramoto [T5]. Its derivation is based on a generalized time-dependent 
Ginzburg-Landau equation. Sivashinsky |32j derived the KS equation as an asymp- 
totic approximation of a diffusional-thermal flame model. The equation also describes 
small-amplitude waves on the surface of a thin film flowing down a planar inclined 
wall (e.g. daHU). 

With the addition of the dispersive term, u xxx , the KS equation becomes the gKS 
equation. Like the KS equation, it has been reported for a wide variety of systems, 
from plasma waves with dispersion due to finite ion banana width [5] to a thin film 
flowing down a planar wall for near-critical conditions (e.g. [29, 18 ). The studies 
in [SJ [3H1 [37] have developed a coherent-structure theory for the interaction of the 
solitary-pulse solutions of the gKS equation. In J6[ [36] the theory was shown to be 
in agreement with experiments using a thin film coating a vertical fiber, another 
hydrodynamic system where the gKS equation can be applicable. 

The well-posedness of (]1.1[) is established for example in [35] in the class of gener- 
alized Burgers equations which consist of a quadratic nonlinearity and arbitrary linear 
parabolic part. The article [20] verifies solvability of the gKS equation in bounded 
domains and studies its limit towards the Korteweg-de Vries equation. In the con- 
text of long-time and large-space considerations, there are recent analytical attempts 
to verify an "equipartition principle" in the power spectrum of periodic solutions by 
deriving bounds on their space average of |u| and certain derivatives of it, see [51 127]. 
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Such a spectral characterization is reminiscent of white noise. 

An interesting work that applies the optimal prediction to the KS equation is 
that of Stinis [33J . Since this approach requires a non- invariant measure, the author 
constructs a Gibbs measure for the required initial distribution through inference from 
empirical data (obtained by a computational approach). This allows then to define 
the conditional expectation providing optimality by an orthogonal projection of the 
unresolved modes to the resolved ones. However, this approach already assumes a 
Gaussian distribution from the outset. For this strategy, one also needs to work with 
the truncated KS equation. Sufficient numerical data are then required in advance 
for a reliable construction of an initial distribution. 

1.2. Notation Functions u£ H s (V a ) for s>l can be represented by their 
Fourier series, 

) = y^feexp( i—x j , u^=u_ fe , (1.11) 



U(X f 

fcez 



where H s denotes here the usual periodic Sobolev space with finite norm, 



\ 2 Hs :=J2(l + \k\ 2 ) S f(k) 



/vt: 



(1.12) 



Furthermore, the square root of the latter quantity is a norm on H s (V a ) equivalent 
to the usual one. We denote for s > 0, 



H s (V a ):=^ueH s (V a ) J udx = 0j. (1.13) 

The subspace of H s (V a ) spanned by the set, 

je^jfceZ, ,-iV<fc<7vj (i.i4) 

is denoted by H%. For a given integer N we define the projections v:=P^u and 
w := Q N u := (I - P N )u by, 



u = PjvM= Ufcexp ( i— x J 

I 7-1 AT \ / 



|fc|<AT 

> = Qnu= u^exp^i— x 



\k\>N 

Let us mention that the gKS equation preserves mass as already noted in Section [TJ 



i.e., 



udx = itQ, (1-16) 



where uq is the zero-th Fourier mode. We remark that Pjv is an orthogonal projection 
with respect to H^j, that means, 

/ (P N u~u)(f>dx = for all <t>€H s N . (1.17) 
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The projection P^r enjoys the following well-known property [131 122] . i.e., for k>s, 
k > it holds, 



\u-V N u\\ H3 <CN 



s — k I 



\ Hk for atiueH%(V a ). 



(1.18) 



Next, we introduce Gevrey spaces. For cr>0 and s>0 we say that a function / 
is in the Gevrey space G CTiS if and only if 



\fi 



^(l + |fc| 2 exp 2^1 + |fcNIM 2 <oc 



(1.19) 



where fk denote the Fourier coefficients of /. Note that if a = 0, then H s = Gq. s . More- 
over, it can be readily be proved, see |17) . that for u<EG atS the following inequality 
holds, 



\u-P N u\\ Hs <N s - k exp(-aN)\\u\\ G 



(1.20) 



2. Formal derivation of a reduced gKS equation As noted in Section [TJ 
we adapt RG approaches [241 25] to the gKS equation. 

2.1. Projections into fast and slow equations We apply the projections Pat 
and Qat defined in (|1.15[) to equation (jl.ip and obtain the following coupled system 
for v and w, 



d t v + P N B(v + w) + A v v = , where A„ =PArA = APjv, 
dtw + QN^(v + w) + A w w = 0, where A w — Q^A = AQjy. 



(2.1) 



We define e = -^, where N is large enough (see error estimates, i.e. Theorem 
and set, 



Ad 



on P N H S =H 



N 5 



A 



jfi on() v //' //' H \. 



(2.2) 



A w - eA v 

The eigenvectors of A„ are the functions exp (i—x) , k £ Z, \k\ < N with eigenvalues, 

k 2 .Jk\ 3 k 4 

Pk '-=-V - -10 - ) +K - . 

a \a J a 

Correspondingly, the eigenvectors of A w are the functions exp (i±x), keZ, \k\>N 
with eigenvalues, 



Pk 



N 4 





k 


2 


'k\*_ 




k 






a 


-iS\ 


K<*J 




a 


1 



(2.3) 



Remark 2.1. The RG method is formally applied here as if these operators were 
independent of e. This technical step of scaling the linear operator and its subsequent 
treatment is part of the abstract RG approach introduced in \24\ WM * n ^ e context of 
fluid dynamics. 

We can now rewrite (|2.1[) by, 



d t v + A v v + P N B (v + w) = , 

d t w + — A w w + QatB (v + w) = . 
e 



(2.4) 
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For convenience we additionally define, 

«-(:)• »-(o). <»> 



and hence rewrite (12 .4[) in the following compact way, 

d t u+-Lu + Au = F(u). (2.6) 

For the subsequent RG analysis we introduce the fast time scale s = | , and we de- 
fine u(s) = u(es). We set u(s) = Pw(s), u)(s) = Qu(s). In this variables (|2.4j) becomes, 

<9 S « + eA„v + eP w B (w + to) = , 
d s w + A w w + cQatB (v + w) = , 

or CT . 

d s w + Lu + e.4u = eF(u). (2.8) 

2.2. Perturbation expansion: The RG equation We now formally apply 
the RG method and additionally omit the dependence of L and A on N as in [24J [25] . 
We make the ansatz of a naive perturbation expansion, 

u e = u° + eu 1 +e 2 u 2 + .... (2.9) 

for u in (|2.8[) . After substituting (|2.9|) into (|2.8|) we formally obtain the following 
sequence of problems, 



(2.10) 



<9 s u + Lm° = 0, 
d s u 1 +Lu 1 =F(u°)-Au°, 

and so on. 

Formally, the solution of (|2.10[h for the initial condition u°(0) = uq is, 

M°(s)=exp(-Ls)w - (2.11) 
Equation (|2.11|) can be equivalently written by, 

v°(s) = v , 
w°(s) = cxp ^-A^sj wo . 

We solve equation (|2. 10|1 ^ with the variation of constants formula, 

u 1 (s) = exp(— Ls) / exp(Ler) [F(exp(— Ler)uo) — -4exp(— Lcr)uo] da, 
Jo 



(2.12) 



(2.13) 



where u 1 (0) = 0, since we are interested in approximations up to 0(e) such that u 1 (0) is 
irrelevant and can be taken to be zero, see [25]. We note that Aexp(— A w <r) —0 (since 
PjvQn =0) such that _Aexp(— Lrr) =A and we decompose the rest of the integrand in 
(I23B1 as, 

exp(Lcr)F(exp(-Ler)u ) -Av =■ F R (u ) + F NR (<r,u ) , (2-14) 
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where Fr(u ) represents the part independent of a on the left hand side of (|2.14l) and 
Fat 7? the rest. Using standard RG terminology, we refer to as the "resonant" and 
Fnr as the "non-resonant" term. 

Using (|2.1ip . (|2.13|) . and (|2.14j) in (|2~9l) provides the following Duhamel's form of 
the formal perturbation expansion for u — u e , 

u e (s)=exp(-Ls) (u Q + esF R (u ) + eJ F NR (a,u ) da^j +0(e 2 ) . (2.15) 

The key idea is now to remove the secular term esF R (uo) which grows in time. To 
this end we define the "renormalized function" U — U(s) as the solution of, 

f-°=*"to. (2 . 16) 

U(Q) = u . 

The equation for the slow variable U(t) = U(t/e) correspondingly satisfies, 

w>-*m, (2 . 17) 

17(0) = n . 

Let us derive the explicit form of the RG equation for our problem. With the expres- 
sions for L and F, and the identity u = v + wq, we get, 

e x p( L.,F ( exp ( -L,) U0 ,^x P (-L.,(:^(» + «^ ( :^; : »). (2.18, 

Next we identify the resonant terms, i.e., F^(u ). With the Fourier series expansion, 



= VWp (i—x) <f> k , 



ke 



(2.19) 



we have, 



B{(j>,ip) = iA^exp ( i^-x ) fc ^exp ( i^-x J 



■ 1/ 



As a consequence, we end up with the expressions, 



(2.20) 



Q N B(v ,exp(-L<r)w ) = i\ ^ exp (i—x j ^ I t; ofe - j exp(-crp ; "')u; o; 



, a / * — ' \ a 

|j|>JV x 7 k+J=j 
/c|<JV<|2 

a 



QjvB(exp(— Lcr)w ,exp(— L<t)w ) —iX ^ expfi — xj (2.21) 



\3\>N 

I 

a 



^ (exp(-crp^)w ofc - j exp(-ap™)w i, 
k+l=j 

\k\,\l\>N 
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\i\ 4 )= 


— K 1 2- 1 







The resonant terms in the first sum are the terms for which pf = p™ holds. This means, 



(2.22) 



lowing set characterizes the resonant indices, 

R 1 (j)~{(k 1 l)\k = Q,j = l,\l\>N} 



The condition p™ + pf = pj characterizes the resonant terms in the second sum of 
(|2.2ip . i.e., (£)" + (^)"= (i;) n for n = 2 ,3,4 needs to hold at the same time. Since 



i I I fc 



" + 2 



k |2 i 



, the set of resonant indices is defined by, 
R 2 (j):={(k,l)\k = 0,l=j,\k\,\l\>N}u{(k,l)\l = 0,k=j,\k\,\l\>N}=$, (2.23) 
since |fc| , |/| > AT. These considerations determine the resonant part of F by, 



-P N B(v )- A v v 
-QatBx(wo,wo) 



(2.24) 



where Bi is given by its Fourier series expansions for the corresponding index set 
Ri(j), i.e., 



Q N B 1 (v ,w ) = 2i\ ^2 c * 

\j\>N 



V 0^ J "u., ■ 



(2.25) 



Equation (|2.14[) and the above consideration give the non-resonant term by, 



F NR (a,u Q ) = 



-P N B ( v + e- A ™ a w ,v + e-~ A ™ a w ) +P N B (v ,v a )+A v v 
-Q N B 1 (v 0l wo)-Q N B 2 (wo,w ) 



, (2.26) 



where Qjv-Bi and QnB 2 are defined by their Fourier series expansions, 



Q N B 1 (v ,w a )=i\ ^2 ° lJaX 

\j\>N 



k+l=3 
|fe|<A r <|i| 



\j\>N k+l=j 
\k\,\l\>N 



J 



a 



vok— )wm + wqi— v k e XPj 



^(pJ-pT) 17 



w k— woie 
a 



{pj-Pk-pY)° 



With (|2.24j) the RG equation for our problem is in the fast time scale, 

d s V + ek v V + eP N B (v,V^j =0, 

d s W + eQ N B 1 (V,W) = 0, 

or after rescaling by t = es, and denoting V — PnU, W — QnU, 

d t V + A v V + P N B(V,V) = 0, 
d t W + Q N B 1 (V,W) = 0. 



(2.27) 



(2.28) 



(2.29) 
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Remark 2.2. 1) The above considerations for the resonant and non-resonant terms 
can easily be extended to general pseudodifferential operators P{A) with symbol 
of the form, 

Rep(iO>c|er, KI-XJO, (2.30) 



where v>3/2. One only needs to adapt the sets for the resonant indices, see \2.22}) 
and V2~2~3\) . 

2) Note that the V -equation in the RG equation V2.29\) is simply the Galerkin approx- 
imation of the gKS equation 

The special structure of the renormalization equation (|2.29jl -? for the unresolved 
(fast) variable allows to give an explicit expression for its solution. After rewriting 
(|2T29]) 9 by 

d t W j {t) + 2i\i-V {t)W j (t)=0, (2.31) 
a 

where Vo(t) — const, due to conservation of mass (|1.16p . we immediately obtain the 
solution, 

W] {t)=c> w e^ V °\ r n : U ; :,.K (2.32) 
With (12321) the solution of (ET2l?|) becomes, 



W(x,t)= 4v e 



i-t(x+2\V t) 



\j\>N 



(2.33) 



Equation (12.33)) shows that there is no restriction on the definition of the mass Vb . In 
the context of stochastic mode reduction the situation is different, see Section 2) 

2.3. Construction of approximate/renormalized solutions In order to 
define renormalized solutions we have first to determine the non-resonant term 
FAffl(c,Mo) given by (|2.26[) . In fact, we are interested in, 



{s,U)= / F NR (a,U)da. (2.34) 



Let 



NR\ 



PF m ( s ,[/) = MVe^ V LlLv k ^Wi 

\j\<N k+l= 3 rt 

\k\<N<\l\ 

\k\,\l\>N 

„ .j „ e {p7-p?) s -\ j ( 2 - 35 ) 
IjI>jv k+i=j p i Pl a 

\k\<N<\l\ 

mm 

e (p™-Pfc -pT) s — i j 



ay*** y — —-Wk—Wi . 

k+l=j 
\k\,\l\>N 



m>N k +l=] Pj~ Pk ~ Pi <* 
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Renormalization of the gKS equation 



Now, we are able to define the approximate solution suggested by the RG theory. We 
obtain, 

u^t)= C -^(U{t) + eY NR (t/e,U{t))), (2.36) 
or with respect to fast w e and slow variables TJ £ , 

= P N u e = V{t) + ePF NR (t/e,U{t)) , 



w < = Q NU e = e-Q« At {W{t) + eQF NR {t/e,U(t))) . 

We note that the initial data are defined by, 

v e (0) = V(0) + ePF NR (0,U(0)) = v + ePF NR (0,u ) , 
w - (0) = W (0) + eQFjv^O, U(0)) = w Q . 



(2.37) 



(2.38) 



3. The renormalized gKS equation and approximation error After in- 
serting (|2.36p into (jl.ip we obtain the following perturbed gKS equation, 

d t u e + Au e + B(u e ,u £ ) = -eR e (U(t)) , (3.1) 
where A and B are defined by (|1.3|) and R c is given by, 

R e ^B^ L W(t)^F NR (t/e,U(t))^+B^e-^F NR (t/e,U(t)),e-^U(t)^ 

+ B(e~ h iF NR (t/e,U(t)),e-^F NR (t/e,U(t))^) - APF NR {t/e, U(t)) (3-2) 
-e- L UuF NR (t/e,U(t))d t U. 



Next, we study estimates on the approximate solutions it e of equation (13. ip . In 
a first step, we need to investigate the non- resonant part Fjvi? of the approximate 
solutions. 

Lemma 3.1. Letp>2 and let g € H q (V a ) with g>4. Assume that the solution of the 
RG equation V2.29]) satisfies U{t) £ H p (V a ) for all t>0. For N large enough there 
exist two constants c\ and C<i, where C% depends on the initial conditions, and c\ 
depends only on V a > but both independent of N , such that the following estimates are 
true for all t>0, 

\\P N F NR (t/e,U(t))\\ HP < C 2 e^ NH , 
\\e-^ M Q N F NR (t/e, U(t))\\ Hp < C 2 e- C - NH . 

Proof, c, c\ , and C 2 represent generic constants independent of N (or e) . We first 
derive estimate (|3.3[h . With the expression (|2.35P i we immediately obtain, 



\\PNF NR (t/e,U(t))\\ HP <ce~ cN * (||V-VW||^ + ||W- VV\\ HP + \\W> VW\ 
where we used the fact that we have the following bound, 



Hp J > 

(3.4) 



PT l/ N * { -„\L\<- i5 (Ly + K \-L\*) 



<ce" ciV \ (3.5) 
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The last inequality follows due to TV < [I]. The second estimate (|3.3p 9 can be obtained 
in the same way by using the inequalities 1 — e~ x < x for all x > and xe~ x < ^ for all 
x> 0. We refer the interested reader to [53] for a deeper consideration. □ 

The bounds of Lemma 13.11 allow us to control R € in the spirit of [51] . 

Lemma 3.2. For N > and g <= H q (V a ) for q > 4, i/iere exist two constants C\ and C% 
independent of N , such that the following estimate holds true for all t>0, 

||i? e WII L2 <C 2 e- Cl7v4t . (3.6) 



Proof. The proof follows in the same way as the proof of Lemma 13.11 We only 
need to take into account the expression of R e and apply Lemma 13.11 □ 

Lemma 3.3. For < T* < oo and g as in Lemma \3.2\ there exists an < e* < oo such 
that for 0<e-~l/N 4 <e* solutions to equation |E2J) satisfy u e (z L 00 (0,T*;L 2 (V a )) H 
L 2 (0,T*;H 2 (V a )). 



Proof. The lemma follows by formally testing equation (|3.1|) with u e and using 
periodicity oiV a , i-e., 4 (9 x (u e ) 3 ,l) =0, such that, 



^IKf-Hl^lf^""' 



i||R,(t/)|| 2 + CM||«'f + «||^«-|| 3 



After defining 



we multiply (|3.7|) by exp ^— / '/3ds^ such that 



Since 



we can rewrite (13.91) as 



(3.7) 



1 2 (3.8) 

7 :=-||i? £ (C/)|| 2 , 



Iexp(-/3t) | || u e || 2 < exp(-/Ji) £ || u £ || 2 +exp(-/3t) 7 (t) . (3.9) 



^exp(-/3t)i|| U l 2 U-/3exp(-/3t)i|| W l 2 +exp(-/30i||| U £ ||, (3.10) 



^ (expH^KH 2 ) <exp(-/3i) 7 (i), (3.11) 

and subsequent integration together with Lemma 13.21 gives. 

hl 2 (T)<Cexp(/3T)^ T exp((/3-C/£)t)^<^^exp(/3T). (3.12) 

For aribtrary < T* < oo we can choose < e < e* := cxp ^T'-)-i3/cc wri ere the con- 
stant Coo is chosen such that, 

^ 7 -exp(/3T*)<C 0O <Cexp(/3T*)//3. (3.13) 
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Renormalization of the gKS equation 



□ 

The reduced equation for the resolved (slow) modes v e alone follows immediately 
after using (jOTlh and (EHU)i, i.e., V(t) = v e (t) - ePF NR (t/e, U(t)), 

d t v e +A v v e +P N B(v e ,v £ ) = ~eG t {U{t),v e ) , (3.14) 

where the induced force term G £ is defined by, 

G e (U(t),v') :=P N B (e- L h e ,e- h <PF NR (t/e,U(t))) 

+ P N B(c-^PF NR (t/e,U(t)),e-^v e ) 

v ' (3.15) 

+ PArB(c- L 7PFAr fl (i/e,[/(t)),e- L ^PFAr fl .(i/ e ,f/(t))) 

-KPF NR (t/e,U(t))-e- h UuPF NR (t/e,U(t))d t U , 

and U(t) = V(t) + W(t) = v e {t) + W(t) is the solution of the RG equation (J2J25). 

Lemma 3.4. For < T* < oo and g as in Lemma \3. 2[ there exists an < e* < oo such 
that for 0<e = -^<e* solutions to {3A$ satisfy v e :=P N u e £ L oo {0,T*;L 2 {V a ))n 
L 2 (0,T*;H 2 (V a )). 

Proof. The proof is similar to the proof of Lemma 13.31 □ 

The same arguments as those for Lemma T3.2I lead to the following. 

Lemma 3.5. For N > and Vq = P]yg € H q (V a ) for q > A, there exist two constants C\ 
and Ci independent of N , such that the following estimate holds true for all t>0, 

\\G t {t)\\ L -<C2^ NH . (3.16) 



Estimate (|3.16j) shows that the RG method provides an invariant manifold by the 
solutions v e such that all orbits are exponentially attracted. The following theorem 
gives qualitative information about the RG approach by quantifying the error between 
(l3~T4l) and (fTTTTt . 

Theorem 3.6. Let g£H 4 {T> a ), e = j^, and suppose that u,u e €L°°(0,T;H 2 (V a )). 
Then, the difference between the reduced solution v e and the exact solution of the gKS 
equation U.l\) satisfies the following error estimate, 

h-v e f LHVa) (T)<Ce 2 +eMGT)(e 1 / i + eJ . (3.17) 

If we suppose that u,u e € L°° (0,T;G (Tj 2{'Pa)) > then we can improve \3.11ty in the fol- 
lowing way, 

||«-» e |ll a(7 , )(T)<C e a + exp(CT)(e 1 /*exp(--5j)+c)). (3.18) 



Remark 3.7. The exponential growth in time is not surprising, see for example 
estimate 13.12\) in the proof of Lemma \3.3i This estimate motivates the definition of 
a new variable h(x,t;r]) := exp(— rjt)u(x,t). Tadmor verifies in f35j/ global existence for 
such a decayed variable h and a conservative form of the KS equation. 
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Proof. The error \\u — v e \\ L2 ^-p \ can be bounded using the triangle inequality by, 

\\ u - v lL'(r a )^\\ u - u lvy> a )Hv 6 - ve \\L'y> a ) > (3-19) 

where the first term on the right-hand side in p,19|) represents the approximation error 
from the RG method (RG error) and the second term accounts for the truncation 
error ( Tr error) . For notational brevity, we introduce the error variables 

e e RG :=u — u e , and e e Tr :=u e — v e . (3.20) 

Step 1: (RG error) The equation for the error variable e RG reads, 

d t e RG + [nd 2 x + 5dl + vdt] e RG + e RG d x u + u*d x e RG = eR e (U) . (3.21) 

First, we test (13.21)) with — d x e e RG , i.e., 

dt (d x e RG ,d x e e RG ) - {nd 2 x e RG ,d 2 x e RG ) - {Sdle e RG ,dle RG ) - (vdte RG ,d 2 x e RG ) 

- (e RG dxU,d 2 e RG ) - (u e d x e RG ,dle RG ) =- (eR € {U),d 2 x e RG ) . 



Then, we use the test function e RG , 

dt{e e RG ,e RG ) + {Kdle RG ,e e RG ) + (Sd x l e e RG ,e e RG ) + (vd%e e RG ,e e RG ) 

+ (e RG d x u, e RG ) + (u e d x e e RG ,e RG ) = + (eR e (U) , e e RG ) . 



(3.22) 



(3.23) 



Next, we add up (|3.22|) and (|3.23[) and apply the Sobolev embedding theorem and 
standard inequalities to end up with, 



1*. 

2~dt 



l4 G ir + ||9xe^ G irj+(^-3a)[||^ G ||^+||^ G | 
<C RG (K,a,e,\\u\\ H1 ,\\d x u\\ H1 ,\\u e \\ H1 ,\\d x u e \\ H1 )\\e RG \\ H - 1 
After defining 

Crg = 2C RG , 

we can multiply (|3.24p by exp (— J Q * Crg ds^j such that 



1 



exp - / C RG ds 



dt 



-RG I 



< exp - / C RG ds 



Crg 



-RG\ 



(3.24) 



(3.25) 



(3.26) 



Applying a corresponding identity based on the product rule as (|3 . 10[) in the proof of 
Lemma l3~3l we can simplify (|3.26j) to, 



<0, 



(3.27) 



which further reduces by assumptions of Theorem 13.61 and after integration to, 



^\\eR G \\ (T)<-exp[ I C RG d s )\\e RG r(0)<C\\e RG r(0). (3.28) 



16 Renormalization of the gKS equation 

In order to get a bound controlled by e on the right-hand side of (|3.28p . we have to 
take into account the definition of the initial data (|2.38l) . i.e., 

fell 2 (0) = e 2 \\PF NR (0,g)\\ 2 < Ce 2 , (3.29) 

since g£H 2 (V a ) and hence PF NR (0,g)eH 1 (V a ) and its norm is bounded indepen- 
dently of e = -jra in 77 1 ('P Q 0- Hence, we can conclude that, 

\\eR G \\ 2 LHVa) (t)<Ce 2 , (3.30) 

which holds uniformly in time. 

Step 2: (TV error) We derive an estimate for the error variable e e Tr := u e — v t — 
u e — PNU e . From (|3.1|) and (|3.14p . the error e t Tr satisfies the equation, 

dte^+lKdl + Sdl + iyd^e^ + PNiu'u'J-v'v^eiPNR^-G^U^)^), (3.31) 

which can be rewritten for all cj) £ Hfj by 

d t (ef, r , 0) + ( [nd 2 x + + ^ f ) 



+ ({Pat - , 4,0 = 6 ({PAri? e (C/) - G £ (C/, , 
Choosing </> = e t Tr allows to estimate (|3.32[) in the following way, 



(3.32) 



^l|eT^| 2 + HI^J 2 < C («^)ll^l| 2 + «ll^,l| 2 + a) + W, ( 3 - 33 ) 



where we define, 



(7) :=(Pjv -OB^ey , 
(J/J^edlP^^II + IIG.^^IDHefvl 



(3.34) 



Let us first control term (7), that means, 

\(I)\<\(-PNU%P N u e ) x + u'u e x ,e e Tr )\ + \{P N u e (P N u e ) x ~v e v e x ,e e Tr )\ 

<^\(d x ((u e -P N u e )(u e + P N u e )),e e Tr )\ + 1/41(9, (P N u e + v e ),e e Tr )\ (3.35) 

<C(h1 H1 + ||PAr W l H1 )ll« £ -P^1 ffl |l4,ll+G||e^|| 2 , 

where we used the embedding H^{V a ) into L°°(P a ), i.e., ||<9 a: (P;vu e + u e )|| i oo < C. 
The second term (77) immediately becomes, 

|(77)| <eexp(-ci/e)C || e e Tr \\ . (3.36) 

Define 

1 {t):=c{\\u^\\M\\^\\ H1 ^\p N u^\ H1 .]\u^\\ H2 )e 1 ^ 

+ C(|| U e ||,||w e ||)eexp(-Ct/e), (3.37) 
(3:=2(C[a,K) + C), 

and mulitply (I3.33P with exp ( — J* (3 ds^j such that 

|exp(~/?i) ^ ||e^ r f < |exp(-/?0 ||e^|| 2 + 7 (t)exp(-/3i) . (3.38) 
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Using again a corresponding identity to (|3.10[) in Lemma 13.31 we can rewrite (|3.38l) 
as, 



d ( , _ 1 



iJt ^exp(-^)^||e^|| 2 J < 7 (t)exp(- / 9t), (3.39) 

which becomes after integration with respect to time, 

-exp(-/3T)||4 r || 2 (T)<C J {(e 1 / 4 + eexp(-Ct/ e ))exp(-/3t)}dt 

< ^ (1 - exp (/3T)) + j^Lj- (1 - exp (- {fi + C/e)T)) . 

(3.40) 

In the remaining part we want to improve (|3.40[) with the help of Gevrey spaces. To 
this end, we remark that the factor e 1 / 4 in (|3.37|) relies on the interpolation estimate 
(|1.18p . If we assume that solutions u e are in G a ,s, then we improve (|1.18|) by (|1.20p . 
As a consequence, we are able to rewrite inequality (|3.40[) by 



i||e^|| 2 (T)<C^ T ( e 1 / 4 exp(/3(r-t)--^)+ e exp(-Ci/ e + /3(T-t))) 



< ^exp (/3T) + ^f^exp ((3T) 



dt 

(3.41) 



□ 



4. Stochastic mode reduction In this section the renormalized equations 
(|2.29|) . (|2.31|) from Sections 12.11 - l3l allow for a rigorous stochastic mode reduction 
similar in spirit to the Mori-Zwanzig one [551 EH] but for systems not satisfying an 
extended or generalized Hamiltonian structure |41) and without an intrinsic invariant 
measure, as is required for Mori-Zwanzig. 

As in the Mori-Zanzig formalism, we assign a stochastic process to the unresolved 
modes. This is done by applying Jaynes maximum entropy principle, see [151 I14[ [28]. 
This seems a reasonable approach for our problem since we do not have a canonically 
induced probability density. Hence maximizing the information entropy for the prob- 
ability density of Fourier modes is equivalent to maximizing the multiplicity of Fourier 
modes. Multiplicity means the number of different ways a certain state in a system 
can be achieved. States in a system with the highest multiplicity can be realized by 
nature in the largest number of ways. Hence, the probability density functions with 
maximum entropy are optimal statistical descriptions. 

It should also be noted that a system at equilibrium will most probably be found 
in the state of highest multiplicity since fluctuations from that state will be usually 
too small to measure. The probability distribution may also be obtained from ex- 
periments as statistical data. In [33], the probability distribution is constructed by a 
conditional expectation obtained from previously computed samples which are used 
to fit an a priori assumed Gaussian distribution. 

Finally, we emphasize that the maximum entropy principle can also be applied 
to problems where one lacks deterministic data as a consequence of not enough ex- 
perimental data to fix all degrees of freedom. A common approach to model such 
uncertainty is to use white noise. The maximum entropy method turns out to be an 
attractive alternative because we do not have to assume the probability distribution. 
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Renormalization of the gKS equation 



However, since we apply the entropy maximization principle |15j on an approx- 
imate equation, we already neglect information from the beginning and hence have 
to account for this by an asymptotic in time characterization of the fast modes for 
example, see Assumption (A) below. This assumption might be improved or adapted 
appropriately in other applications. 

4.1. Problem induced probability density by maximizing information 
entropy With the considerations at the beginning of Section |4j we assign a proba- 
bility distribution to the unresolved degrees of freedom W based on the following 

Assumptions: 

(A) For the unresolved degrees of freedom W we define the weighted variance 
E[V(W)]by, 



E[V(W)}:=E 



12 

|fc|>JV 



(4.1) 



for arbitrary cr^ > and (3 £ R. We assume that the fast modes already 
asymptotically reached their equilibrium state in time, that means, 



E[8 t V(W)]=ai t p - L , (4.2) 

for alH > large enough. 

(B) Under (i) the probability density function P[Wjv+i <Wn+i,Wn+2 < 
w N+2 ,...] :=F(wm+i,wn + 2,...), i-e., 



F(w N+1 ,w N+2 , ...):= / / ...f(w N+1 ,w N+2 ,...)dw N+1 dw N+ 



has maximum information entropy Si(f), 

5 J (/) = -//(w)]og(^)dw > 



2 • • • , 

(4.3) 
(4.4) 



where dw := dwN+idwN+2 ■ ■ ■ , /( w ) denotes the probability density of the 
unresolved variable W and i/(w) is an invariant measure which is defined on 
background information intrinsically given by the physical origin of W. 

Remark 4.1. 1) The idea of deriving probability distributions for multiscale evolution 
problems by maximizing the information entropy seems to go back to \21f . The energy 
argument in \21}/ . which assumes that the fast modes reached already the stationary 
state, does not provide here enough information to fix the Lagrange multiplier Ai 
associated with this energy constraint. We impose Assumption (A) instead. 

2) A mechanical system governed by the Hamiltonian H(q,p) canonically induces an 
invariant measure by the density distribution function f(q,p) :— zW) e ~^ H p ^ ■ 

3) Equation \4->fy accounts for the fact that the fast modes do not entirely contribute to 
an invariant measure. For simplicity, we also neglect a possible random contribution 
to the evolution in time. This is a further reason for the strong assumption in fi4-2\ ). 

In information theory, an entropy related to (|4.4[) was originally introduced by 
Shannon [30] to measure the maximum information content in a message. The As- 
sumptions (A) and (B) above account for the lack of a free energy and a Hamiltonian 
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for which the thermodynamic equilibrium (invariant measure) can be achieved via the 
gradient flow with respect to the Wasserstein distance [T6|. In fact, it should be noted 
that minimizing the free energy with respect to constant internal energy is equivalent 
to maximizing the entropy. 

Theorem 4.1. Under Assumptions (A) and (B), it follows that the unresolved modes 
Wk for \k\ > N obtained by equation 12.31]) are normally distributed with zero mean, 



-ff 



i.e., /ife = 0, and variance Ofc = 2 \±v > where Xi(t) := is a Langrange multiplier. 

To keep the considerations simple, we only account for a spatial random process 
and keep the time deterministic in Theorem 14. II (and Assumption (A)). 

Proof. To maximize (j4.4[) under Assumptions (A) and (B), we apply the following 
constraints: 

(CI) l//( w ) dw = 1 ' (45) 
V ' {E[d t V(W)}:=Jf(w)d t V(w)dw = a 2 00 ^-\ V " ' 

where (|4.5[> 9 is a consequence of assumption (i). 

Since a Hamiltonian structure for our problem does not exist as already noted, 
it is not clear how to define the invariant measure v under the constraints (|4.5|) . We 
will see in equations (14.14[) and (14.17[) that the following simple choice for v will be 
sufficient, 

-( w ) : = II (4.6) 

\k\>N k 

Remark 4.2. The measure v{w) is the probability density function if we only have 
a priori information. Usually, it is a nontrivial task and basic considerations of 
symmetries are required to find this measure v. 



Hence, maximizing the entropy Si subject to the constraints (|4.5I) leads to 

£/(w){log( J + f(w) + Ao + AiflkV(w))dw = 0, (4.7) 



where A^ for £ = 0,1 are Lagrange multipliers associated with the constraints (|4.5[) . In 
order to give (|4.5|> o a precise meaning, we write down the explicit form of the equation 
belonging to each Fourier coefficient of the fast mode variable w — Qn u solving (|2.29|) . 
We briefly show what the constraint (I4.5[> 9 means, 



E[ftV(W)] = ^//(w) 



WkdtWkdw 



\k\>N 

12 v r (4 - 8) 

fc|>JV 

Remark 4.3. If we suppose zero mass, then we obtain for the variable Vq representing 
the zeroth Fourier mode in 12. 3 1]) . (2.32]) , and (2. 33]) that Vq =0. As a consequence, 
the solution of the renormalization equation of the unresolved (fast) variable W(x,t) 
is defined by its initial condition in this case, 



W(x)=w (x), 



(4.9) 



20 Renormalization of the gKS equation 

see 12.31]) . This raises the question of whether it might be not more suitable to gen- 
erally work in the case of mass Vq = 1? 

Since (|4.7p should hold for arbitrary variations Sf, we obtain the following ex- 
pression for the probability density function, 

/(w) = li/(w)e- Al * v(w) ) (4.10) 

where Z := c I/ ( w )+ Ao is called the "partition function" which is determined by the 
normalization constraint (14.51) 9. i.e., 



Z:=Z(X 1 )= / ^(w)e- AiatV(w) dw. (4.11) 



Since the constraint (14.51) 9 is quadratic in its nature, we represent it by 

12 ^ XiVo , 12 1 1 



w(w)= _ ^ _ wfe= _ ^ ___ (u , fe _ Mfc ) , (412) 

where 



|fc|>JV '" " " |fc|>iV" W 



aw = olTTF ' and A^ = 0. (4.13) 

With identities (|4.12[) and (|4. 1 3[) the probability density function (|4.10[) can be written 
as 

\k\>N 

where M is the normal distribution given by 



^11 c Wk ka w V ZirN (p k ,k<j w ,w k ), (4.14) 



7V(/i fc) fca w ,^) = - ^=e , (4.15) 

fcuwv 2vr 

which is characterized by the following moments 

POO 

Af((j,,a,w)dw = l, / Af(n 7 o-,w)wdw = /i, and 



(4.16) 



A/"(/i, cr,w)w 2 dw — a 2 + fi 2 



The first property in (I4.16|) together with the normalization condition ()4.5j) i allow us 
to define the partition function Z by 

Z= II <w k ka w V2n- (4,17) 

|k|>JV 

The probability density function / admits then the simple form as a product of Gaus- 



sian distributions, i.e., 



/( w )= N{v.k,kaw,w k ). 



(4.18) 

|fc|>JV 
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With the second and third property in (|4.16p and the constraint (|4.5|) 9. i.e., (|4.8p . we 
obtain for all \k\> N, 



The important information contained in formula (|4.20[) and (|4.13|) i is that we do not 
have to assert to each Fourier mode k its standard deviation fc 2 <r|y. We only need 
to determine once the Lagrange paramter Ai via (|4.20[) . We already know from the 
beginning ( (|4.19[1 ) how to choose the mean, i.e., =0 for all \k\ >N. □ 

Hence, the approach of maximizing the generalized information entropy allows 
to systematically determine the probability distribution function /(w) of the Fourier 
modes for the unresolved degrees of freedom W. 

We emphasize that the RG approach suggests multiplicative noise as a com- 
pensation for the unresolved modes unlike the commonly obtained additive noise by 
Mori-Zwanzig's mode reduction [4"T] , Moreover, an estimate (13.61) . which can be cor- 
respondingly derived by additionally accounting for the Galerkin error, shows that 
the influence of the stochastic force decreases for decreasing e := -A^. 

5. Direct approach: Replacement of G e by white noise The result 
of Lemma 13.51 also enables for a direct approach to model the unresolved degrees 
of freedom as completely unknown. Such kind of complete uncertainty is generally 
described by white noise W(x) with zero mean and a variance equal to the power 
spectral density. It is very common and widely accepted to model uncertainty by 
white noise. Hence, we replace eG e (U(t),v e ) in equation (|3.14p by 



where W(x) & L 2 (V a ) is the Gaussian random variable as motivated above, i.e., with 
zero mean /i and suitable variance a. It is immediately clear that M e is a compatible 
replacement of G e since (|5.1[) satisfies a bound corresponding to the one in Lemma 
13.51 One can follow Stinis' approach [3 3) for example in order to determine fi and a 
by a maximum likelihood method. 

6. Discussion and conclusions We have formally developed a new stochastic 
mode reduction strategy with a rigorous basis by obtaining appropriate error esti- 
mates. The analysis can be summarized in three key steps as follows: 

(1) RG method: The RG technique [231 [25] turns out to be a formal and feasible 
method to decompose the gKS equation into slow v e and fast variables w e , respectively. 
The equation for the slow modes v e represents a Galerkin approximation of the gKS 
equation plus an additional perturbed force term eG e which also depends on the 
infinite dimensional renormalized fast modes W. An important property of the RG 
technique is that it can be easily extended to higher space dimensions, see [24]. We 
also remark that the dispersion term, i.e. u xxx , does not affect the mode reduction 
analysis. 




(4.19) 



We conclude with (|4.19[) that the Lagrange parameter Ai is, 




(4.20) 



Af e (x,t) := eexp(-Ctye) W(x) , 



(5-1) 
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(2) Error bounds: We rigorously characterize the formal RG method (1) by qual- 
itative error estimates fTheorem l3.6l) . These estimates further allow for an additional 
direct mode reduction strategy which is much simpler and straightforward but not 
as systematic. The basic idea is to replace the perturbed force term eG e directly by 
white noise. A physical motivation for such a simplified reduction is the fact that 
white noise is a well-accepted random model for complete uncertainty. 

(3) Maximum entropy principle: Due to the lack of a Hamiltonian structure and 
an invariant measure, we apply Jaynes' maximum entropy principle [151 114] to define 
the renormalized fast modes W as a random variable. This random variable then, 
together with the renormalized approximation of the slow variable v e , provides a 
systematic explanation for the appearance of a noisy low dimensional gKS equation. 
In contrast to optimal prediction we obtain optimality in the sense of maximum 
entropy here. 

There are three main features of the new low dimensional gKS equations: 

(i) Reliable and efficient numerics: The low dimensional formulation developed 
here should allow for reliable (since information from the unresolved degrees of free- 
dom included) and efficient (since low dimensional) numerical approximations. In 
fact, we systematically account for the unresolved degrees of freedom by the steps 
(1) and (2) above. This is especially of importance since the choice of slow and fast 
variables depends on the physical problem and is often not clear. For instance, by 
considering the gKS in large domains, it is possible to introduce a further scale which 
accounts for the unstable modes. Hence, we can study three different scales such as 
"unstable modes" , "slow stable modes" , and "fast stable modes" . The main question 
is then how to account for the unstable and the fast stable modes in an equation for 
the resolved slow modes only. 

Moreover, the error estimates from step (2) provide a qualitative measure on how 
to choose the dimension of the slow variable. This is also the main advantage of mode 
reduction considerations over pure convergence analyses of Galerkin approximations 
(e.g. numerical schemes) where one completely neglects the unresolved degrees of 
freedom. Hence, straightforward discretization strategies might loose model relevant 
information in the neglected degrees of freedom. This is a major motivation to in- 
clude rigorous mode reduction strategies as an important part of the development 
of computational schemes. We further remark that this is a major reason why the 
application of noise to deterministic partial differential equations shows good results 
and is currently a topic of increasing interest. 

(ii) No Hamiltonian structure; no invariant measure: Many classical mode re- 
duction strategies rely either on a Hamiltonian structure or an invariant measure. 
Based on the three steps (l)-(3) above, the new asymptotic reduction strategy cir- 
cumvents such dependences. For example, when classical optimal prediction methods 
[5] fail because of such deficiencies, the stochastic renormalization provides optimal- 
ity in the sense of maximum information entropy and hence proves as a promising 
alternative. 

(Hi) The role of noise: We gain a rigorous understanding of the origin of noise 
and the way it appears in the gKS equation. This is especially of interest due to 
numerical evidence provided together with a heuristic motivation in [4] for instance. 

Clearly, there are open questions and future perspectives. For example, motivated 
by the comparative study initiated by Stinis |34j , it would be of interest to numerically 
analyze and compare available mode reduction strategies such as adiabatic elimination 
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[38) . invariant manifolds [TJ, and optimal prediction [3] with the new RG approach 
developed here. Since the statistically based optimal prediction [3^ is performed for 
a truncated KS equation, it provides a convenient setup for comparison with the new 
and more generally applicable method suggested here. 

Another question is how can we apply the RG method to the derivation of a low 
dimensional approximation for a gKS equation investigated under three scales, i.e., 
"slow unstable modes" , "slow stable modes" , and "fast stable modes" . 

The RG method is based on a natural splitting into linear and nonlinear terms by 
the variation of constant formula. Recent studies, e.g. by Holden et al. [HUE], make 
use of such a splitting via a suitable numerical scheme for equations with Burgers' 
nonlinearity. Hence, the reliability and efficiency of the renormalized low dimensional 
gKS equation motivate the application such numerical splitting strategies to the new 
reduced equations derived here. Finally, we emphasize that efficient low dimensional 
approximations are of great interest for numerical scrutiny of long time asymptotes. 
We shall examine these and related issues in future studies. 
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